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ABSTRACT 


This thesis presents an analysis of a vivaldi antenna. It 
is a broad band antenna realised in a printed form. The most 
commonly used numerical technique. "The moment method" is 
discussed. Its implementation is explained by applying it to 
some simple structures like a dipole and a rectangular microstrip 
antenna. An IBM-PC compatible A-D card is fabricated and used to 
measure the approximate form of current distribution on the 
antenna. The theory of the sampling probe is presented. Finally 
the antenna is made to fit in an elliptic cylinder coordinate 
system and approximate form of bases functions given based on the 
experimental results- 
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CHAPTER I 

INTRODUCTION 


INTRODUCTION 

The developments in the field of antennas have come a long 
way through from the conventional antennas like dish antenna, 
yagiuda antenna etc. to the modern microstrip antennas. With the 
advent of high speed computers and new mathematical tools to 
solve the difficult electromagnetic equations numerically, the 
antenna field has taken a great stride. The advantages of 
microstrip antennas being numerous like light weight, low cost, 
easy reproducibility etc. CIU all the applications of 
conventional antennas are being tried to be handled by the patch 
antennas. One such innovative step is the realisation of broad 
band antennas in printed form- All the microstrip planar 
antennas are of the form which has a resonant cavity and hence a 
low bandwidth, but the introduction of vivaldi antenna by Gibson 
et.al. C2I3 has led to the use of printed antennas in broadband 
appl i cations . 

LITERATURE SURVEYi 

Until the introduction of Vivaldi antenna by Gibson CZJ and 
end fire tftpered slot antenna by liahapatra and Prasad C33 there 
was no well known planar antenna which could produce a well 
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defined end fire beam with appreciable bandwidth. Gibson C23 
showed that the vivaldi antenna is capable of producing symmetric 
beam over a wide band width His antenna consisted of 
metallized dielectric substrate with an exuponent ial Hot in the 
metallization. This type of antenna is capable of multi-active 
bandwidth- The length and operative size of such an antenna are 
determined by the lowest frequency of operation. At the lowest 
frequency length has to be greater than wavelength, and the 
operative size of about half wavelength for efficient radiation. 

R-K. Janaswamy C43 has analysed the antenna to calculate the 
radiation pattern by assuming it to be approxximated by a stepped 
slot line. The geometry of the antenna and the approximation are 
as show* in the Fig. C S- 3. 

But this method was valid only for antenna whose operative 
size is small. His method of analysis consisted of two steps. 
In the first step, the tangential component of the electric-field 
distribution in the tapered slot, is obtained. In the second 
step, far fields radiated by the equivalent magnetic current in 
the slot are obtained using an appropriate Green’s function. 

The aperture distribution in the tapered slot is determined 
by employing the usual travelling wave antenna assumption that 
the aperture distribution is governed predominantly by the 
propogating modes. He divides the antenna into a number of 
sections of line of uniform width connected end to end. Then use 
is made of the theory of small reflections to get an overall 


reflection coefficient. 




F,5 C. Geometries as solved 
in [4] and Is'] 
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Later the same author C5D solved the antenna by 
approximating it as shown in Fig. C i- 3. 

The solution was based on the Moment Method Solution of the 
electric field Integral equation for the electric surface 
currents flowing on the conducting plates of the antenna. The 
geometry is such that use can be made of simple rectangular 
patches in the Moment Method Solution. The presence of the 
dielectric is included in the Moment Method Solution via the 
impedance boundary conditions. The current distribution and 
Input Impedance were calculated. The method is valid for 
arbitrary aperture widths. 

S.N. Prasad and Mahapatra C3I3 gave the experimental results 
of a slot line antenna that operates in X-band with a gain of 6 
dB and a side lobe level of 10 dB. 

Ehid Gazit C63 gives a novel method of feeding the vivaldi 
antenna. He explains about the effects of transition from 
miscrostrip line to asymmetric slot line, and the effect it has 
on the band width characteristics. 

The various attempts to solve the vivaldi antenna had moment 
method formulation using Subdomain functions. The use of 
Subdomain functions has the disadvantage of taking a lot of 
computer time. Even the use of surface patch dipoles as 
subdomain functions took a time of 8 to 10 hours C*r 3. Also the 
shape of the vivaldi antenna is such that ordinary subdomain 
functions like rectangular or triangular patches cannot be used 
to represent the antenna completely. The approximations which 
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have to be made in the shape of the antenna are many due to which 
the analysed antenna cannot truly represent the actual vivaldi 
antenna. 

AIM OF THE PRESENT WORK: 

Keeping in view the difficulties involved in the analysis of 
the vivaldi antenna, an effort is made to find the reasons for 
the drawbacks involved, the complexities of the moment method and 
an attempt to solve it in a different fashion- 

In Chapter II the design details of the antenna, the feed of 
the antenna and other particulars discussed. The formulation of 
the Integral equation and the most commonly used numerical 
technique, the moment method are also discussed. The moment 
method is applied to a dipole and a rectangular patch microstrip 
antenna to discuss about the complexities of the method. 

In Chapter III the significance of current distribution and 
the experimental determination of the current distribution on a 
vivaldi antenna are discussed. Also the insertion loss 
characteristics of the antenna are given. 

In Chapter IV an attempt is made to confine the antenna in 
an orthogonal coordinate system and better way of dividing the 
antenna into subdomains discussed. Approximate form of the 
subdomain functions are discussed. 

In Chapter V summary and scope for further work in this area 


are discussed. 
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CHAPTER II 


MOMENT METHOD 


INTRODUCTION: 

The vivaldi antenna consists of an exponential asymmetric 
slot line- CFig.S-H- The antenna is etched on Cu clad Duroid 
sheet. The antenna radiation is in the end fire direction with 
the E— field in the x direction. The slot width of the radiating 
end and the length of the tapered slot line section are kept 
approximately same. When frequency is such that the slot width 
is approximately half wavelength, the radiator has high 

effeciency and can work effeciently beyond this frequency. The 

transition is made from a microstrip line to an assymetric slot 

line over a length of around half a wavelength and hence the 

input reflection coefficient increases as frequency is decreased 
below- The Fig. C 2- H also defines the coordinate axis for 
further analysis. 

For a symmetric slot line CFig.^ □ the fields are given by 





l-ig. 2- Geometry of the fabricated vivaldi antenna 







E = -r? H 

' X r 


where 



€ +1 
r 


is the wavelength in the dielectric 


and rj = Intrinsic Impedance of air. 

For jx| large (:<) and (j<) may be written as 

• I I % e~l 

H (j :< ) = 

" rr^ 

when I X j = r . 

It can be seen from the above equations that the higher the 
dielectric constant, the better will the slot line confine the 
fields very near it and function as a good transmission line. 
Hence a low dielectric constant material has to be chosen to make 
the tapered slot line work effectively as an antenna. The value 

. tt 

of the dielectric chosen is 2.Z and the thickness is 

The advantage of the assymetric slot line antenna when 
compared to the symmetric slot line antenna which is present only 
on one side of the dielectric is that the input transition is 
better beause of the low impedance of the assymetric slot line, 
which can be easily matched to a microstrip line. 

Also the double sided slot antenna can be used more 
effectively while designing a circular array for an 


omnidirectional antenna- 


The two major subjects of vivaldi antenna design ares 
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(a) The transition from the main transmission line (usually 

microstrip) to a slot line, for feeding the antenna; this 
should have a wide operating frequency range and low 

reflection coefficient. 

(b) The dimensions and shape of the antenna, to obtain the 

required beamwidth, bandwidth, sidelobes and backlobes. 

The microstrip, the double parallel strip and the slotlines 
are all of 50 Q impedance. A vivaldi antenna is fabricated and 
is used as a termination to the transition. 

In order to evaluate the radiation pattern and the near 

field, the exact current distribution is required. Also the 
Input Impedance as a function of frequency can be evaluated once 
the current distribution is known. In order to find the current 
distribution the boundary conditions are satisfied on the antenna 
surface and an Integral equation obtained, which is later salved 
by a numerical technique such as the moment method. 

FORMULATION OF THE INTEGRAL EOUATION: 

The key to the solution of any antenna or scattering problem 
is a knowledge of the physical or equivalent current density 
distribution on the volume or surface of the antenna or 
scatterer. Once these are known then the radiated or scattered 
fields can be found using the standard radiation integrals. A 
main objective then of any solution method is to be able to 
calculate accurately the current densities over the antenna. 
This can be accomplished by the Integral Equation (IE) method. 

In general there are many forms of Integral equations. Two 
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of the most popular for time harmonic electromagnetics are the 
Electric field Integral equation (EFIE) and the magnetic field 
integral equation (fflFIE). The EFIE enforces the boundary 
condition on the tangential component of the electric field, 
whereas the NFIE enforces the boundary condition on the 
tangential components of the magnetic field. 

The EFIE is based on the boundary condition that the total 
tangential electric field on a perfectly conducting surface of an 
antenna or scatterer is zero. This can be expressed ass 
= E^^ + E^® =0 on s 

or E^^ = -E^^ on S. 

where S is the conducting surface of the antenna or scatterer. 

E^ = Total tangential electric field. 

E^^ = Incident tangential electric field and 

E^^ = Scattered tangential electric field. 

The incident field that impinges on the surface S of the 
antenna or scatterer induces on it an electric current density J 

s 

which inturn radiates the scattered field. If J is known, the 

s 

scattered field every where that is due to can be found using 
the relation 


E^(r) 


- j 


/J -e 


V (V . A) 


(A ^ 




w p«A+V(^.A) 


where 


A( r ) 


h JJ "s 




dS’ 


hence 
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where , 


+ V JJ J^(r*) G(r, rM 



G(r, r’) 


4fT ^ 


^-j/? \r - r^\ 
4 n 1 r - r’ j 


R = I r - r’ j 

If the observations are restricted on the 
antenna (r = r^ ) then the above equation becomes 


surface 


of the 


2 [ J J (r) G<r^, r-) ds- 

s 

^ JJ ^ 3^<rM G(r^, 7’) ds ’ j = (r = r^) 

Because the right side of the equation is expressed in terms 
of the known incident electric field it is referred to as the 
electric field integral equation (EFIE) . It can be used to find 
the current density 3^ (r') at any point r = r' on the antenna or 
scatterer. Once the integral equation is obtained it is solved 
numerically by a method such as the moment method. 


THE MOMENT METHODi 

There are various methods to obtain the exact current 
distribution, prominent among them being spectral domain 
technique, immittance approach, transmission line modelling, 
moment method etc. All of the above methods and many more can be 
used to obtain the solution of the integral equation which is the 
basic equation having the information about the current 
distribution . 

Analytical solution of the integral equation is rarely 
possible to obtain, and in almost all the cases numerical 



techniques are used, the most common being the moment method, 
C103 knowing the voltage at the feed terminals and determining 
the current distribution, the input impedance and the radiation 
pattern can then be obtained. 

An integral equation is of the form 


F(g) = h — (t) 

where F is a known linear (integral) operator, h is a known 
excitation function and g is the resonse function. The objective 
here is to determine g once F and h are specified. Though the 
solution is often not possible in closed form, the linearity of 
the operator F makes a numerical solution possible. The 
technique known as the moment method requires that the unknown 
response function is expanded as a linear combination of N terms 
and is written as 


g(z') C^g^ (zM + (z') + ... C^g^ (z’) = 

Each C is an unknown constant and each g (z*) 
n n 

function referred to as basis or expansion function. 
( 1 ) reduces to 


N 

r C g (z* 
n^n 

n=1 

is a known 
Equat ion 


2 c„ F(g^) = h _W 

n=1 

The basis functions g are chosen so that each F (g ) can be 

■^n ^n 

evaluated conveniently, preferably in closed form or at the’ very 

least numerically- The only task then is to find the C ; the N 

n 

unknown constants. 


BASIS FUNCTIONSi 


One very important step in any numerical solution is the 
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choice of basis functions. In general one choses as basis 
functions, the set which has the ability to accurately represent 
and resemble the anticipated unknown function, while minimising 
the computational effort required to employ it- These basis 
functions are divided into two general classes the subdomain 
function, which are non zero only over a part of the antenna 
surface and the continuous entire domain functions which are non 
zero over the entire antenna. Some examples of the subdomain 
functions being piecewise constant, piecewise linear, piecewise 
sinusoid etc. as in Fig.^ . A common entire domain basis set is 
that of a sinusoidal function. The main advantage of entire 
domain bases lies in problems where the unknown function is 
assumed apriori to follow a known pattern. In such cases, entire 
domain functions may render an acceptable representation of the 
unknown while using far fewer terms in the expansion than would 
be necessary for subdomain bases. 

Eqn. (H) is a linear equation with N unknowns and hence to 
obtain the N unknowns we define the inner product 

< w, g > ~ J J '^-9 

S 

where S is the surface of the structure being analysed and w is 
the weighting function. The equation then reduces to a matrix 



equation 
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where 



The matrix equation may be solved by inversion and can be written 
as 


-1 



As special cases if weighting functions are chosen to be same 
as basis functions, then it is called Galerkin's method and if 
they are chosen as delta functions then it is called point 
matching. 

MOMENT METHOD APPLIED TO A DIPOLE: 

To illustrate the solution outlined above we display the 
results for a dipole of length 0-5 X and diameter = 0.005 X. • The 
subdomain bases functions used are pulse functions and weighting 
functions used are delta functions. We would have assumed the 
current distribution on a dipole to be sinusoidal, but though the 
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magnitude of the current does have a sinusoidal form, the 
standing wave pattern is not completely valid, in the sense that 
there are phase variations in the current along the antenna 
length. Thus the information about the reactive portion of the 
current is obtained by applying moment method. 

The current distribution and the Input Impedance calculated 
for a dipole in Fig.C^'l is shown in the Fig.C^U and Table 

From the Fig. c6d it is clear that the current distribution 
on a dipole is having a real^and an imaginary part. The real 
part is greater in magnitude than the imaginary part and is more 
or less sinusoidal. The real part is the major contributor for 
the far field pattern and this is the reason for obtaining a good 
idea about the radiation pattern of the dipole even if we assume 
a sinusoidal current distribution with a constant phase. Though 
the actual current is not truly sinusoidal with the same phase. 
The sinusoidal distribution with same phase completely fails to 
give zero electric field on the dipole surface because, the 
imaginary part has a considerable contribution towards the near 
field. When the actual current distribution is used then the 
tangential electric field is found to be negligible on the whole 
of the antenna when compared to the field at the 
point. It is found to have the following forms 


excitation 





The values of the input impedance are calculated for two 
different type of source modelling: (a) the deltagap modelling 
and (b) the magnetic frill modelling. 

For finding the radiation pattern accurately the source 
modelling might not play an important role because one deals with 
far field which is not a very sensitive function of current 
distribution, but for the calculation of the input impedance, 
source modelling is a significant factor. The way we model the 
source should be an accurate representation of the actual source. 
For the case of the cylindrical dipole a coaxial fed monopole 
would be better represented at its excitation point by a magnetic 
frill whereas a delta gap would be a better representation for a 
twin wire feed. 

Also since point— matching is a method whose solutions 
satisfy the electromagnetic boundary conditions (vanishing 
tangential electric fields on the surface of an electric 
conductor) only at discrete points, the boundary conditions may 
not be satisfied at the in between points and the deviation is 
known as the residual C103. For a half wavelength dipole, a 
typical residual is as shown in Fig.C-^U. The residual decreases 



if Galerkins method instead of point matching is used. The 
results for Galerkins method are shown in Fig.CT^"!. 

MOMENT METHOD AS APPLIED TO MICROSTRIP ANTENNASi 

A microstrip antenna consists of a metallic plate or patch 
on a dielectric slab CFig.9 1 . The patch is fed either by a 

microstrip transmission line contacting an edge of patch or by a 
coaxial probe extending through the ground plane and contacting 
the patch. 

The microstrip patch when fed by a microstrip line and a 
coaxial proble will have different source models. In the former 
case the microstrip feed line can be replaced by an equivalent 
current source which is obtainable from the transverse electric 
or magnetic fields (E^ or H*) in the plane where the microstrip 
feed line connects to the patch as shown in Fig. C ^ The 

transverse fields can be represented by the TEM feed line fields 

.V _ ^ I 

E, = - z — and = - y 

where w^ is the effective width w^ = w + A w where w is the strip 

width and Aw is an extension which accounts for the fringing 

fields at the edges of the feed line. From uniqueness concepts 
backed by a perfect magnetic conductor is needed; therefore 
the microstrip feed line can be replaced by an electric current 
source y where y = z / -J w^z ’ . Likewise for a probe fed 

patch, the probe can be replaced by an equivalent current source 
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FIG-9 GEOMETRY OF A MICROSTRIP ANTENNA 
AND EQUIVALENT OF A MICROSTRIP 
ANTENNA EXCITED BY A MICROSTRIP 
FEED LINE 
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y present on a cylindrical surface with diameter equal to 

d where w = z / (rr d ) and d is the diameter of the probe, 
o o o 

Moment method can be used to solve the microstrip antenna 
C113. The method is based upon a boundary condition or integral 
equation formulation, with the unknowns being the currents on 
microstrip patches and the wire feed lines plus their images in 
the ground plane. The integral equation is solved using the 
method of moments. The integral equation is formulated with the 
zero reaction concept C123 developed by Rumsey. 

THE REACTION TECHNIOUE: 

Consider the exterior scattering problem illustrated in the 
Fig.CiUJU. In the presence of a dielectric or conducting body, the 
impressed electric and magnetic currents (J^, M^ ) generate the 
electric and magnetic field intensities (E, H). Let the exterior 
medium be free space. 

From the surface equivalence theorem of Schelkunoff C iO 3, 

the interior field will vanish (without disturbing the enterior 

field) if we introduce the following surface current densities 

J = n X H 

s 

M = E X n 

s 

on the closed surface S of the scatterer. In this situation, 

illustrated in Fig.Cdl? 3 we may replace the scatterer with free 
space without disturbing the field anywhere. 

By definition, the incident field ( E^ , ) is generated by 

< J^, M^ ) in free space, and the scattered field is 




N / 


FIG. 10 ZERO REACTION THEOREM 
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E = E - E. 

s 1 

H = H - H. 

s 1 


when the surface current 

< 

M 

5 

) radiates in free space. 

it 

generates the field ( 

E , 

H ) 

in 

the exterior and (— E., H. ) 

in 


s 

s 


1 1 


the interior region. 

The 

result 

is shown in Fig . fd.C>3 

We 


place an electric test source in this region and find from the 
reciprocity theorem that 

n ' - "s ni < Ji Et - Ht Wv = 0 

where is the free space field of the test source. The 

above equation is the integral equation for the scattering 
problem and this can be converted into a matrix equation using 
moment method. We get 


czj cij = cvn 

Typical terms of z and V are as given below: 
Z 


mn 


J E*" . J ds 


n 


n 


V 


J e"* . J . dl 


m 


n 


where the surface current has been expanded as 
N _ 

J = E T J and 
s n n 

n=1 

is the impressed current filament. 


MODIFICATION FOR DIELECTRIC SLABi 


If there is a dielectric slab present in the antenna then 
the slab is removed and replaced by free space and the equivalent 
volume polarization currents 
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~3yj = jw E 

The matrix equation then gets modified to 
£ 1 + A 1 1 CI3 = CV3 

where 

AZ = -jw(fi-€)r dv 

mn 0 J d 

V 

where Ej is the electric field radiated by 3 the n— th expansion 
d n 

mode in the presence of the dielectric slab, and e”* is still the 
free space field of the m— th test mode. For a microstrip antenna 
with ground plane E^ is given in a simplified manner in C13I1. 

For our problem use of Sommerfeld integrals is required for the 
evaluation of E^. The feed model is a very important factor 
determining the accuracy of the results obtained. 

□nee current distribution is known then the input impedance 
and radiation pattern can be easily found. 

RECTANGULAR PATCH MICROSTRIP ANTENNA: 

To illustrate the method outlined above a rectangular patch 
antenna without dielectric has been analysed and the results 
presented. The geometry of the antenna is shown in Fig.JCdJj 
.along with the results. There are serious limitations if 
rectangular pulse bases are used and the maximum cell width is 
restricted to 0.1 X when the size of the antenna is greater than 
2 or more square wavelengths in area. 

Here use is made of sinusoidal reaction formulation as given 
by Richmond C14I1. In this method the surface of the antenna is 
divided into cells and the surface current distribution is 
expanded in overlapping sinusoidal bases in the direction of 
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current flow and rectangular pulse bases in the transverse 
direction as in Fig. Ci2ll via an application of Galerkins method, 
the integral equation formulated with the zero reaction concept 
is reduced to a matrix equation. 



where the subscript b refers to the bottom plate and t refers to 
the top plate. From symmetry of the problem we have 


^b ^tt ^bb' ^tb 


Z . * and 
bt 


-V^ as a consequence of which 


tt tb 


The results for an antenna which is 15 cms by 7.5 cms is 
given in the Fig. CJ33. 

Also the basis function of sinusoidal form CFig 443 over a 
patch has the advantage that it has the exact expressions for the 
near field and hence resulting in reduction of one order of 
numerical integration C14I1 in calculating the impedance elements 
of the impedance matrix. In general the sinusoidal basis 
functions closely approximate many wave forms and also helps in 
reducing the number of elements or order of impedance matrix 
because of their greater domain. 

The rectangular microstrip antenna with a ground plane is 
found to be resonant at a frequency of 933 MHz and a small 


bandwidth. 
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i 

5 . 1 883 

933 . 4 

60 . 8033 

*!■ 

j 

5.2177 

60 . 7538 

■f 

j 

5.8147 

933 . 6 

60.8415 

Hh 

i 

5.2514 

60 . 7920 

+ 

i 

5 . 8478 

933.8 

60 . 8796 

+ 

j 

5 . 2851 

60.8301 


j 

5 . 2797 


Observat ions 

Ztt « self impedance “ Rs + j Xs 
ztb * mutual impedance =* Rm + J Xm 

1 ) Always Rs > RM 

but Rb Rm 

2) Xs > Xm if frequency > 933.0 MHz 
Xs < Xm if frequency < 933.0 MHz 

Xs ■ Xm Somewhere near frequency » 933 . 0 MHz 
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The dipole and the microstrip patch were easily solved 
because of the fact that they were simple in structure and 
confined to a coordinate system. The solutions can be more 
accurately and quickly obtained if the approximate current 
di s t r i bhut i on which is known a priori is used as entire domain 
functions. In the case of rectangular patch microstrip antenna 
use was made of sinusoidal surface patch dipoles which are entire 
domain functions for this problem. 



CHAPTER III 


MEASUREMENT OF CURRENT DISTRIBUTION ON VIVALDI ANTENNA 

INTRODUCTIONS 

Hidden in the integral equation is the information about the 
current and once it is known, it can be used to find the radiated 
fields using the equation 

E^(r) = - jwA-j — V (V .A) 

W >£ 

where A is the vector potential. 

The most important properties of an antenna, the radiation 
pattern and the input impedance are dependant solely on the 
current distribution. Though there can be an intuitive feel of 
the form of current distribution (like assuming a sinusoidal 
distribution over a dipole), the exact value of the current at 
each point is required to obtain the radiation pattern and the 
input impedance. Though an approximate form may give a fairly 
good picture of the radiation pattern, it completely fails to 
give a correct value for the input impedance and the near field 

L73. The knowledge of the exact current distribution is felt all 

/ 

the more necessary when one has to deal with antenna arrays where 
antennas are placed very near to each other and the near field of 
the neighbouring antennas is felt on them. Analytical solutions 
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of the integral equation are possible rarely and can be found 
only for two or three cases of antennas like spheroidal antennas 
where the antenna confines itself to a particular coordinate 
system. If the antenna is of any other shape there is no other 
way than to resort to approximations and numerical techniques 
like moment method to solve for the current distribution. 

The moment method formulation requires the appropriate 
selection of the basis functions and weighing functions. A wrong 
selection of the basis functions and weighing functions may not 
give accurate results unless the number of iterations are many 
which in turn means a lot of computer time. The weighing 
functions play an important role in the sense that they are the 
ones which decide the forcing of the boundary conditions. Using 
weighing functions like delta functions might lead to very poor 
results because the boundary conditions are satisfied at only 
very few discrete points, whereas galerkins method of choosing 
the same basis functions as weighing functions gives better 
results because of the fact that the boundary conditions are 
satisfied on an average over the surface of the antenna and not 
at discrete points. Also the number of terms required will be 
very less for convergence if the entire domain functions are used 
and better results will be obtained if they closely approximate 
the actual current distribution. For some antennas like dipole 
it is easy to choose the entire domain functions as ^sinusoidal 
functions because it closely approximates the actual current 
distribution on the antenna. 
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Hence an attempt is made to get an approximate form of 

current distribution on the antenna by measurement. 

EXPERIMENTAL MEASUREMENT OF CURRENT DISTRIBUTION; 

Since it has been observed that the current distribution is 
the most important fator which determines the radiation pattern 
and input impedance, we will explain methods to find the current 
distribution experimentally. Measurement of the current on the 
surface of an antenna involves measurement of the tangential 
magnetic field very near this surface. These measurements are 

made in the reactive near field of the antenna, and the 

requirement that they be made in very close proximity to the 

antenna may dictate the use of very small probe antennas to 
sample these fields- From an experimental point of view, the 
problem is to design or obtain probes or receiving antennas with 
associated apparatus for determining the amplitude, phase and 
polarization of the electromagnet! c field at any point in space, 
in the immediate vicinity, including the surface, of radiating 
antenna . 

sampling PROBE - THEORY AND DESIGN; 

On a perfectly conducting surface, the surface current 
density 7 is related to the magnetic field by the relation 

n X H = J 

where n is the unit normal directed outward from the surface. 
The depth of penetration of the tangential magnetic field in a 
good conductor is so small in and above the VHF range of 
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frequencies that the true current distribution is very closely 
approximated by the surface distribution on a perfect conductor. 
Thus if the output of a probe is proportional only to the 
magnitude of the magnetic field it may be considered to be a 
measurement of the true surface current. Suitable probe for the 
measurement of the tangential magnetic field is an infinitesimal 
magnetic dipole which to a good approximation may be realised fay 
a small loop antenna. 

The following requirements are to be met by the probe if an 
accurate measure of the current is required CSH. 

1. Probe should not distort the field. 

2. Probe should be so small as to represent a point. 

3. Probe should have the desired polarization 

4. The output should be of a magnitude, capable of being 
measured . 

Simple loop — The loop antenna consists of a single turn of wire 
as shown in Fig. ciS H. 

If the loop is small and the magnetic field uniform then the 
current If, induced in the loop fay the magnetic field is equal in 
magnitude at all points on the contour. The total current 
through the load is 


I, = + I + I 
L H z y 

where is due to the magnetic field and are due to the 
electric field. Due to symmetry does not contribute to the 
current but if present contributes to the load current. Hence 
it is required to keep the loop in such a way that E field is 
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tangential to the symmetric portion of the loop. 

The magnetic sensitivity of small loops is proportional to 
the loop area and hence the loop must be of sufficient size to 
provide a useful output. Balanced against this is the fact that 
the output of the probe is proportional to the average value of 
the magnetic field through the loop aperture, hence the aperture 
must be small enough to essentially measure the field at a point. 
(Also the presence of non-vanishing longitudinal electric field 
at points away from the antenna surface means that the magnetic 
field away from the surface is no longer exactly proportional to 
the current). 

To overcome the above mentioned difficulty double loaded 
loops CFig are used. 

The total current in due to the fields across the two 
gaps can be resolved into symmetric and antisymmetric components- 
The symmetric components I are proportional to and 

antisymmetric components are proportional to E_^ . 

In a practical system Ig and can be obtained by using a 
180 ° hybrid. 






"a " 


where v is the isolation of the hybrid. 


DESIGN OF A-D CARD* 

The measurement of the tangential magnetic field involves 
movement of the probe all over the antenna in a predetermined 
fashion and then noting down the values of the position of the 
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probe, values of amplitude and phase of the received power 
But this recording involves a lot of strain and inaccuracy in 
measurement due to the large number of values to be recorded. 
Also the noted values have to be later entered on to a P.C. to 
get an idea of the approximate form of current distribution. 

The above difficulty has been overcome by designing a PC 
card which with the associated software helps us to record the 
values directly on to the P.C. The total set-up is shown in Fig. 

( ) • The input to the PC Card is from the rear panel of the 

network analyser which gives a voltage of 10 mv/dB for amplitude 
and 50 mv/degree for phase. Thus by properly calibrating the A— D 
converter we can directly get the values of amplitude and phase 
on to the PC which can be stared in a file and later graphs of it 
can be viewed- 

The block diagram of the PC Card is shown in the Fig- 5. 

It can be extended to include more channels so as to get a direct 
measure of frequency also. 

DESIGN DETAILS: 

An extension of PC*s 8038 /jp— bus called I/o channel is 
available on eight, 62 pin edge connectors C9I] CFig. I. Any non 
system hardware, with proper decoding scheme can be put on this 
bus and accessed by the CPU. 

The I/o channel contains an 8 bit bidirectional data bus, 20 
address lines, 6 levels of interrupt, control lines for memory 
and I/o read or write, clock and timing signals, a channel check 
line, 3 channels of DMA control lines, memory refresh timing 



4-5 
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control signals, and power and ground lines. All signals on I/o 
channel are buffered to provide sufficient low-power Schottky 
(LS) loads per slot. 

The 8088 ptp is capable of addressing 64 K different ports 
via its address lines AO — A15. However a PC system board 
decades only lower 10 address lines (AO — A9> restricting the 
number of available ports to 1024. The lower half of the 1024 
ports are reserved for the system board itself, and the upper 
half is dedicated to I/o channel. Many of the I/o channel ports 
are however reserved for the IBM adapter cards. I/o address 
space 300 - 31 FH allotted to IBM prototype card, are used for 
present application. 


A9 AS A7 A6 A5 A4 A3 A2 A1 AO 

110 0 0 

always fixed bits used for decoding 

The inclusion of AEN signal for the address decoding is note 
worthy. This ensures that decoding is not done during DMA 
cycles. The various chips used are 
74LS138 - Decoder 

74LS245 - Buffer 


4052 


Multiplexer 


7574 


A-D converter 


Tfie circuit diagram is shown in Figure 


casi 


MEASUKEMEMT: 

The setup used to measure the current distribution is as 
shown in Fig. The input to the antenna is given by port 1 
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and the received power at port 2 is a direct measure of the 
tangential magnetic field and hence to the current on the antenna 
with some proportionality constant. 

^ I tan H I a Current density 

The probe is held normal to the antenna and is made to move 
along x— axis for different values of y and along y— axis for 
different values of x. The orientation of the probe is possible 
in 2 directions and each direction gives a different component of 
current. Thus when we have the knowledge of the form of the 
tangential magnetic field, we have completely characterized the 
form of the current distribution on the antenna. Since we have 
information about amplitude and phase of the current 
distribution, it tells us whether we have a standing wave or 
travelling wave and also about the nature of decay present. 

The data when observed for different values of x for 
constant y and different values of y for constant x gives us the 
form of current distribution at that cross section. The graph of 
phase for a standing wave pattern would be as in Fig.CjrO 3 and 
for travelling wave as in Fig.C^Ol. Similarly depending on the 
decay the graph of the amplitude part takes its shape, some of 
the shapes being given in Fig. 

RESULTS AND DISCUSSION: 

Measurements have been made at regular intervals of x for 
constant value of y and for regular intervals of y for constant 
value of X- The graphs of amplitude and phase variation at 
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different cross sections are given in the Figs. to Li^l and 

the results discussed in the following pages. 

The total current distribution J can be split as follows: 

J = J(!<, y) = xJ^(x, y)+yJ^{x, y) 

= xJ (x).J (y)+yJ<x).J(y) 

X X y y 

* 

In the following pages approximate mathematical forms are 

given for J (x), J^(y), J (x) and J (y). 

X X y y 
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From the plots of the phase at the 6 cross sections it can 
be seen that the phase has a linear variation with x and hence 
can be said to have a travelling wave distribution. 

From the plots of the amplitude at the 6 cross sections it 
can be seen that the amplitude decreases as the distance of the 
observation point is away from the inner curves (1) and (E) . 
Thus we can assume that J (x) is a function of x with an 

X 

attenuation constant and a pohase constant 








<3L h cdLi^ 4- 



From the plots of phase at the 6 rross section it can be 
seen that the phase pattern is that of a standing wave pattern- 

The amplitude falls rapidly as we go from the inner edges of 
the antenna (1) and (2) towards the outer edges (3) and (4). 

The current distribution J (y) can be said to have the forms 

X 

(a y) 

J (y) = e sin (.f'i y) 

where because amplitude increases as y increases on the 


antenna surface. 
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From the plote of phase at the different cross sections 
can be observed that J^iy) has a standing wave patterns. 

The amplitude decreases as we go from inner curves of 
antenna (1> and <E) towards the outer curves of the antenna 
and (4). Hence the current distribution J^Cy) can be said 


have a form 


(ot y ) 

Jyty) = e sin y) 

a > 0 « 

yy 


it 

the 

C3) 

to 


where 
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The graph of insertion Loss vs. Freq. is given in Fig. €.%S 
D.It can be seen that the insertion loss is less than -10 dB for 
frequencies right from 1 GHz to 12 GHz. Insertion loss of -10 dB 
implies that 90/1 of the incident power is radiated and only 1071 
is reflected due to impedance mismatch. The antenna is resonant 
at a frequency of 7.4 GHz when it has an insertion loss of around 
-40 dB, which means that almost all the incident power is 
radiated. 

If dB insertion loss is considered to represent a good 
antenna then the vivaldi antenna can be seen to have a bandwidth 
of GHz. 

In order to improve the bandwidth and also understand the 
mechanism of radiation etc-, knowledge of current distribution is 


essential . 
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CHAPTER IV 


VIVALDI ANTENNA IN ELLIPTIC CYLINDER COORDINATES 


INTRODUCTIONi 

It has been observed in the previous chapters that the 
vivaldi antenna is of a shape, which cannot be solved directly by 
moment method using rectangular or triangular subdomain basis 
functions and gross approximations had to be made to solve the 
antenna. To get results nearer to that of an actual vivaldi 
antenna, we will try to confine the antenna in a coordinate 
system which closely approximates the shape of vivaldi antenna in 
two dimensions. Also since we have got a fairly good idea of the 
form of current distribution from our measurements, we will use 
them to get better basis functions. 

ELLIPTIC-CYLINDER COORDINATES [151: 

It is obtained by taking an orthogonal family of confocal 
ellipses and hyperbolas in a plane and translating in the 
z-direction Fig.CS^D. The coordinate surfaces are elliptic, 
cylinders (ry = constt), hyperbolic cylinders (y/ = constt) and 
parallel planes (z = constt). The relation between the new 
coordinates and rectangular coordinates may be written as 




Fio. Klliptin-ryliiulor coordinnlon. Flin cdordinato ;;iirnu'c'5< nro olliplic 

cylinclcrB y) “ liypcrbolic cyluidera «// = const^ and parallol planes z = const. 
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K = a cosh n cos V' 

y = a sinh r/ sin v 

z = z 


The surfaces n - conslt are elliptic cylinders. 

2 2 
X y 

2 + 2 = 

(a cos h r/ ) (a sin h 77) 

while the surfaces v' = constt. are the hyperbolic cylinders. 



ENTIRE DOMAIN FUNCTIONS* 

Till now the methods used for solving an antenna for its 
current distribution always preferred subdomain basis functions, 
the reason being that the coefficients can then be obtained to 
match any unknown current distribution, the antenna has- But 
witl* the advent of arrays in wliich thousands of elements are 
present and especially where the antennas are very close to each 
other, the time taken to account for all the array element 
currents and also the mutual impedance effects in addition to the 
dielectric character izat i on needs a lot of computation time. 
Also the use of subdomain basis functions can be better only for 
atitennas whicli are small in terms of wave 1 engtl^s . 

But once the form of current distribution is known, then use 
can be made of entire domain functions, whicft offer various 
advantages in addition to very less computation time. 



6y^ 

MeasuremGnls wfiiich were made on vivaldi antenna gave the 
current d i s t r i !m t i on whirls was oKplaiiind in tl^e previous chapter . 
From tf)e graphs it can be observed tl^at the currer^t distribution 
followed a particular pattern whf»n we try to analyse the vivaldi 
antenna in tl>e elliptic cylinder coordinate system^ we need to 
have tlie l)asis functions in tl^e two orthogonal directions a^ and 
a 

We can see f r ()rn tfie diagram l)elnw that tive ry current wfiich 
enters at the input tends to more along the curved path of the 
aritenna and a rougl) approximation would be to assume that for - 
constt, the current has a travelling wave distribution with a 
decay which accounts for tl>e radiation field. 

It can also be seen that 77 current has standing wave 
distribution for r/ conslt berause it can he ap|>r cnc ima ted as a 
wave front of equal phase. 
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Our experimental results also show a similar behaviour. The 
above analysis then allow us to assume entire domain functions 
which have the mathematical form: 

J = r/ (Y7, vw) + V (y), 

Y) J (Y7 ) . J (V ) + \w J (y>). J (vw) 

rt y) ^ 

where (>7), <V')f (t) > y (V^) has the forms as was just 

eKplained, with attenuation constant accounting for the decay and 
phase constant accounting for the phase transition- 

In order to solve for the current distribution on the 

vivaldi antenna, we divide our antenna into a number of 
subsections as shown in Fig- Zjjl-lr on which we assume only y) 

directed current distribution for simplicity and assume that on 
each strip it is constant in the v' direction. 

For finding the impedance elements we need to have 
expressions for near field in the elliptic-cylinder coordinate 
systems . 

EXPRESSION FOR NEAR-FI ELD-EVALUATIONt 

If we assume that only y) directed current is present and it 
is of the form 
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1 

J V X {1 

R J 

>7 ^ < 0 • . 

V' ’ , ^ ‘*'7 ’ 


V X 

(gF) = 

(Vg) ;< F ~ g (V k 

F) 

H = 

1 

4fT 

Jv 

R J 

XT) f (o’ 

, V' ' » z ' > dyj » 

= 

1 

4rT 

1 { 1 

l+j/^R . 

r2 j 

e-j/7R J 

« ■[ R 1* il7 ’ , yf,- ^ 2 . ) J tty) ' 

The expression for relation between unit vectors in 

elliptic-cylinder coordinates and cylinder coordinates is given 

in Appendi x 1 ) . 

Fi = 

1 


(x - X* ) 

+ j (y-y* 

) C j 


IH-JflR 

-j/^R ! 

r 

[ » 9- ail, V'- 


1 

J 

L / 

'3n> 


+ j 

sinh y) * 

cos v'jj 

f {4« 

(y) r 2 cjy) • 


( 9 „) 


(x-x') sinh y> * cos ' 


-k (y-y’) cosh r/ ’ sin v' * - i (z-z’) sinh yj • 




1+j/^R 


j (z-z*) cosh r) * sin v" 


• ] ] 


f # iyy’r V', z’)} <*7- 


It has 4 components. Consider only 1 component i 


E = 


j <«j « 


V X R 
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1 


j 

1 


V X J g F dr,’ 


j -<= 

1 ♦ j/?R 


V X J - i (z- 2 ’> sinh r, ’ cos v^’ 


(gii> 


1/2 


-j/'j'R 


[ — j f 4 # <ry • , , z ’ > dr) • 


Again using the identity 


Vx(gF) = (7g)xF-g(^xF) 

we can expand the above expression- Similarly for the other 3 
componerits . 

Finally on simplification we get an expression of the form 


E = 


f (a U + a U„ + a 
J r, 1 H z 


U ) dn* 


where 


U = f I y) ’ , z’» r), y, z 

Since we now have a single integration formula for calculating E 
field we can get the impedance elements 

Z = - r E . J ds 

mn J m n 

and thus we cajt fill tl*e impedance matrix- We can tf,en find the 
current distribution by inverting the matrix and multiplying it 
by excitation vector- 

CI3 = CZ3"'’ CV3 

where CI3 = unknown vector 

nZ^ = impedance matrix 


CV3 


excitation matrix 
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As a convenient check it can be used to find the current 
distribution at the input point where it should have a form 
similar to that of a assymetric slot line 

Since the antenna is solved in its natural coordinate system 
we get a better feel of the current distribution and also of tl»e 
radiation mechanism. Also the computation time will be decreased 
because of the use of entire domain functions, because of which 
various antennas with different n V'g limits can be 
evaluated for their bandwidth and other properties. This idea 
can be extended to the numerous antennas which are confined to 
the different coordinate systems like using sinusoidal patch 
dipoles for rectangular antennas. 
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CHAP! ER V 

SUMMARY, CONCLUSIONS AND SCOPE FOR FURTHER WORK 

SUMMARY AND CONCLUSIONSj 

In this thesis an exponentially tapered slot antenna known 
as vivaldi antenna is analysed. The antenna is fabricated on a 
low dielectric substrate. This antenna is different from other 
microstrip printed antennas. The usually present ground plane 
and a resonant cavity between the top surface and the ground 
plane are absent. The antenna is a tapered transmission line 
which radiates as the slot width increases. The taper is 
extended over a lengtl^ of X/2, thus realizing a large bandwidth. 
Also the transition from microstrip line to the antenna is made 
over a length of X/2. The antennas insertion loss 

characteristics were measured and found to show a large 
bandwidth . 

The most widely used moment method to solve the integral 
equation is discussed and applied to a dipole antenna and a 
microstrip patch antenna. The advantages of entire domain 

functions and tl»e significance of the source modelling are 
discussed. The knowledge of exact current distribution is found 
to be necessary for finding out the near field and also the 
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radiation patterr* and input impedance. Input impedance is found 
to be sensitive to tHe way we model tbe source. 

In order to obtain a good idea about the form of current 
distribution, measurements of the near field are done- The forms 
obtained are discussed. Finally the antenna is confined to an 
elliptic cylinder coordinate system and is divided into 
subdomains on which current can be assumed with the help of the 
measured patterns. Since ttie antenna is solved in a coordinate 
system in which it confines to some orthogonal surfaces it can 
lead to better results. 

SCOPE FOR FURTHER WORK: 

Further theoretical work needed are the programming of the 
proposed method for various valves of n’s and v'*s and the 
evaluation of the bandwidth for these antennas. Study of the 
transition from microstrip to assymetric slot line and 
minimisation of the length of tl»p transition are also important. 

Mutual couplirjg behaviour has to be formulated to study the 
performance of elements in array environment. This is very 
important because at spacings small compared to wavelength, the 
mutual coupling is very strong and will determine the input 
impedance behaviour with frequency. 

Finally tl>e design of a circular array of such anterrnas with 
a radial power divides C 16 I and combines to get tl>e desired 
omnidirectional power pattern can be carried out. 
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APPENDIX 1 

RELATION BETWEEN THE UNIT VECTORS OF ELLIPTIC - CYLINDER 
AND RECTANGULAR COORDINATES 


The conversion of a vector from one coordinate system to 
another coordinate system requires the knowledge of the 
relationship of the unit vectors of the coordinate systems. The 
unit vector i is perpendicular to the surface n = constt and 
similarity n is perpendicular to » = constt. The combining 
relation between the E coordinate systems is 


X 


a cosh rj cos 


y = a sinh y) sin v' 

a 


V X = 


i (x) = i 

dx 


V (a cosh n cos y) 


r/ 


a 


(gi,) 


1/2 a-r, 


(a cosh y> cos y) 


where 


'11 


'22 




(g,,> 


a 


1/E a>t.> 


(a cosh n cos 


hence 


X = 


V 


cgii) 


s 


1/2 


(a coshy) cosy-') 


a 




a 


, ,1/E a^ 

( q. . ) 


(a coshy) cosv'l 


'11 
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Similarly, 


(a cosh r/ cos 


(a sinh n cos y' ) 




‘See’ 


Solving equations (1) and (2) we get. 


i = '■ — Y/z sinh n cos \y 1 cosh rj sin y/) 


— — <a cosh r> sin y ) sinh >7 cos \t')J 
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APPENDIX 2 

SOFTWARE FOR DATA ACQUISITION 


Program atd; 
Uses dos, crt; 
Const 


ad-amp = *0318 
ad-phase = *0319 
ad— rd-low = *031A 
ad-rd-high = *0318 
ad-cs = *03 1C 
zero = *00 
ESc = # 27 


Type 

DataArray = array CO.. 5013 of byte; 
VoltArray = array CO.. 5013 of longint; 

Var 

outfile ! text; 

hour, min, sec, se clOO* wor d ; 

por t— Addr , I , J , channe 1 : intege r ; 

dummy - byte ; 

ch ! char ; 

Amp, phase sDataar ray 
Vol t s Vo 1 tar ray 

Function Sample -ATD ( port-Addr : integer > : byte; 
Var 


dummy : byte 



begin 


ff 


dummy := port Cport-AddrD 
dummy != port Cad-rd-low3 
dummy5= port Cad-rd— highU 
delay ( 1 > 

dummy := port Cad-rd-Iow3 
delay ( 1 ) 

sample-ATD:= port Cad-csH 

end; 

Begin 


Writeln ( 'A to D' ) 
writeln; writeln; 
ch:= ’ • ; l:= 0; 

geltime (hour, min, sec, sedOO); 

Writeln ('time:’ haur:3, min::3, sec: 3, 5ec100:3); 

Writeln ('lo start press key*); 
ch : =r eadkey ; 

Repeat 


l:=I+1 

Write 

( ’Voltage=’ ) ; 

Volt(i):=i; 

AmpC13 :=sample-ATD (Ad— Amp); 
phase CI3:= sample-ATD (Ad-phase); 

Writeln (*Amp!=', AmpC13, ’phase:*, phase C13:6); 
Writeln; 

Writeln (’press any key (ESc to end)’; 

Writeln; 
ch : =r eadkey ; 

Until ch=ESc; 

assign (outfile, ’Data out’); 
rewrite (outfile); 

Writeln (outfile, ’x-axis (mm) phase Amplitude*); 

f or j : = 1 to i do 

Writeln (outfile, voltCJ3:6’, ’ ’, AmpCJ3: 6, phae CJ3:6); 

Writeln (outfile ’Go to lotus’) 
close (outf i le ) 
end. 



